Lattice hyper-solitons in photorrefractive materials 
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We show a novel kind of nonlinear waves in two-dimensional photonic lattices. This waves take 
the form of light clusters that may fill an arbitrary number of lattice sites. We have demonstrated 
by numerical simulations that stable propagation can be achieved under adequate conditions and 
we have described the unstable patterns developed otherwise. Our results show that these new kind 
of nonlinear waves can be easily found in current experiments. 
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PACS numbers: 

Optical lattices are artificial crystals of light P] that 
can be obtained by shining an adequate optical material 
(usually photorrefractive) with a number of mutually co- 
herent laser beams. The interference pattern created by 
those beams, induces a spatially periodic modulation of 
the refractive index in the material. This effect can be 
used to trap a probe laser beam, usually much weaker 
and of different frequency than those beams used to gen- 
erate the lattice. 

The study of nonlinear waves in periodic potentials has 
become an active field of research in recent years due to 
its applications in control and manipulation of light prop- 
erties. In this kind of systems, the appropriate combina- 
tion of the probe beam power and the lattice depth and 
period yields to stable and localized light propagation for 
the probe laser. This beam, under adequate conditions, 
can generate a lattice soliton^, linked to a site of the 
periodic potential, that propagates undistorted. 

The generation of lattice solitons has been demon- 
strated experimentally in both one-dimensional 0,0 and 
two-dimensional 5, 6] potentials. Novel kinds of nonlin- 
ear waves like soliton dipoles and quadrupoles0, as 
well as discrete vortices have been also demonstrated]^ 
[Tol |. Most of this results have been generalized for the 
case of several probe beams that are mutually incoher- 
ent, yielding to the so-called discrete vector solitons 11]. 

In this Letter, we will show that all these previous 
distributions belong to a more general and novel kind of 
nonlinear waves that exists in two-dimensional periodic 
lattices. This waves are clusters of parallel laser beams 
that are linked to the sites of the lattice. The adequate 
choice of beam powers and lattice parameters yields to 
the formation of arrays of lattice solitons that propagate 
without shape distortion for arbitrarily long distances. 
The extension of this light distributions can be arbitrary 
large, filling an unlimited number of lattice sites, so they 
may be regarded as hyper-solitons. 

We will calculate the explicit form of these distri- 
butions and demonstrate that stable propagation is 
achieved for a wide range of configurations. We will also 
test the stability of these structures when they are gen- 
erated using arrays of Gaussian beams IT3l| with ade- 
quate powers, widths and relative phases, calling in that 
way for the experimental demonstration. 

We will study the paraxial propagation along z of a 



laser beam with envelope 'i'{x,y,z), through a photor- 
refractive material with a periodic modulation of the re- 
fractive index in the transversal plane {x,y). This pro cess 
can be described by the following wave equation 
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where is the transverse Laplacian operator, the spa- 
tial variables x and y are expressed in units of the wave- 
length (A) and z is measured in units of 47rA. The func- 
tion V{x,y) = Acos^{-Kx/T)cos^{-Ky/T) defines the peri- 
odic optically-induced lattice, where A and T are respec- 
tively its amplitude and period. Both parameters can 
be easily controlled by acting on the laser beams that 
generate the grid. 

We look for stationary solutions of the form 
'^{x,y,z) — u{x,y) ex.p{il3z), where /3 is the propaga- 
tion constant. These solutions can be numerically cal- 
culated and constitute families of solitons described by 
the parameter f3. In Fig. ^ we show some examples of 
stationary states for the particular values of the parame- 
ters, ^ = 1, T = 10. In order to have localized solutions 
with u{x, y —^ ±oo) = 0, due to the particular form of 
the potential term in Eq. it is required that (3 < 0. 

Each family of solitons is topologically defined by the 
number of lobes and phase configuration. In Fig. 12 we 
show plots of the power P = \u\'^dxdy versus the 
propagation constant for several families of solitons, cor- 
responding to different cases illustrated in Fig. ^ 

In the linear limit, the periodicity of the medium re- 
quires Bloch-typc periodic solutions (see Fig.QjA), which 
extend over the whole lattice. If the nonlinear effect is 
significant, localized solutions appear, taking the form of 
multi-hump solitons whose lobes are located at the max- 
ima of the lattice. It can be appreciated that families of 
solitons with a higher number of lobes present a higher 
power. For a particular family, the shape and amplitude 
of the lobes depend on the value of the propagation con- 
stant, so that P basically increases with /3. 

In Fig.n}B,C it is shown that as power raises, two close 
lobes of the same sign undergo fusion into a single lobe 
state (Fig. ^JD) centered at a lattice nodal line, which 
may reach unlimited high values. This is due to the non- 
linear effect which shades the lattice for high powers. 
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FIG. 1: Plots of the amplitude of different multi-hump soli- 
tons, illustrating the dependence of their shapes with power 
(or /3). Dashed lines indicate the nodes of the lattice. Labels 
correspond to points in Fig. |5| 
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family, increasing jS yields to a fusion of the same-sign 
doublets (Fig. ^G) ending with the state in Fig. ^F. 

Fusion and splitting of lobes always correspond to an 
inversion in the sign of the slope of the P vs /? curve 
and this inverted-slope section of the curve is more pro- 
nounced as many lobes are involved in the fusion (see 
for instance thick and thin continuous lines in Fig|21-(b), 
corresponding to two and four lobe family respectively). 
Analogous fusion or splitting take place for states with 
diagonally-disposed lobes (see Fig.^I). In this case, the 
fusion takes place towards a fundamental soliton cen- 
tered at a node in the common corner of both lattice 
sites (Figs.[nJ,K). 

In cases where an odd number of symmetrically dis- 
posed lobes does not permit a preferential fusion of two 
of them, a new scenario appears. For example the family 
of soliton M (Fig.O, changes when (3 is increased, low- 
ering the amplitude of the central lobe up to the point 
where it disappears. Soliton N, however, does not quali- 
tatively change. 
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FIG. 3: Examples of higher-order solitons, showing different 
configurations, including one with a large number of humps 
(labeled R). For all of them P — —0.5. Dashed lines indicate 
the zeros of the lattice. 



FIG. 2: Power curves for different families of multi-hump soli- 
tons, (a) Solitons with two side by side light spots (continuous 
thick line: two same-phase lobes; dashed line: same but in- 
verted phase; continuous thin line: two double-hump spots 
with inverted phase), (b) Solitons with lobes of the same 
relative phase and disposed in different configurations (thick 
dashed: one lobe; thick continuous: two side by side lobes; 
thin dashed: two diagonally disposed lobes; thin continuous: 
four lobes disposed in a two by two array). Letters label 
points corresponding to the examples in Fig. Q Point D, out 
of the frame, has coordinates (—0.02, 1.45 x 10^). 



Modes with two inverted-phase lobes (Fig.QjE) exhibit 
a different behavior with power. Neighbor lobes prevent 
each other from growing up, increasing power arbitrarily 
with (3. Instead, the lobes drift away from each other, 
also increasing their width, up to the point where they are 
located on lattice nodal hues (Fig.QjF). There, the power 
curve reaches an end point where branches corresponding 
to other families join (see Fig.[21-(a)). This is the case for 
the family of solitons composed by four lobes disposed 
in a two double-lobe configuration (Fig. Qjll). For this 



States with a higher number of lobes present more com- 
plicated dependences of the shape with power, and they 
can be basically accounted in terms of the same effects 
described above. An example is the two-by-two lobe soli- 
ton (L) whose lobes also merge together at the negative 
slope region of the curve. On the other hand, branches 
corresponding to different solitons may join together at 
points where fusion of lobes produce the same final soli- 
ton. This is the case at point (K) where the two and four 
lobe solitons fuse into a fundamental one. 

The stability of those hyper-solitons was investigated 
by means of numerical simulations using a standard 
Beam Propagation Method. The evolution of the maxi- 
mum amplitude for some states is plotted in Fig. ^ The 
first conclusion is that those states whose lobes alternate 
sign are completely stable (for instance, all those corre- 
sponding to the dashed line in fig. [2Ia)). Lines E, P, Q 
and R in Fig.0J-(a), correspondent to the same label soli- 
tons in Figs.^andlSl confirm the stability of those states. 
The numerical simulations were also carried out using ar- 
rays of Gaussian beams instead of the exact eigenstates 
and they resulted stable as well. This robustness would 



3 



eventually allow to reproduce the results experimentally 
and make this beams suitable to be used for optical con- 
trol operations. 
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FIG. 4: Maximum amplitude versus propagation constant for 
a number of states. Line labels correspond to the solitons in 
Figs. [Hand 01 

Contrariwise, solitons with lobes of the same sign or 
different sign pattern develop a number of instability 
scenarios which can be summarized as follows. For low 
power, the beam lobes are centered at lattice sites and 
they maintain their position developing power oscilla- 
tions due to the coupling taking place from one to each 
other (see line B in Fig.^(b); solitons H, L and M present 
a similar behavior) . For states corresponding to the neg- 
ative slope region of the power curve, close lobes merge 
together and oscillate laterally with an amplitude of the 
order of the lattice pitch. This transversal oscillations 
also affect the amplitude (line C in Fig.^(b); similar be- 
havior is found for soliton G) . Similar oscillations are also 
found for solitons in the region of the power curve where 
the slope turns again positive, for moderate power. For 



high power, however, the lobes are found to fly off from 
each other after some propagation distance. 

A different scenario is that of states with lobes dis- 
posed diagonally. Alternating sign ones remain stable. 
For those with same sign lobes, however, the instability 
pattern is different according to the branch of the power 
curve. Those belonging to one of the branches (solitons 
1 and N for instance), develop a periodic phase (sign) 
inversion between alternate lobes (see line I in Fig. ^ 
(c)). Nevertheless, solitons of higher order develop for 
longer propagation distance more complicated coupling 
patterns (see line N in Fig. 2|-(d)). On the other hand, 
states belonging to the other branch (like soliton J) de- 
velop a fusion of lobes followed by transversal oscillations, 
similarly to the case of negative slope. 

Concerning the fundamental soliton (one lobe) it is sta- 
ble for low powers when centered at a lattice maximum 
(see simulations for different powers in Fig. QKe)). For 
higher powers, however, it wanders around, since the lat- 
tice is not enough to clamp it to the site, flying away 
at the end. The propagation distance at which the lobe 
starts to move increases as power increases and the soli- 
ton becomes virtually stable as power tends to infinite. A 
similar behavior is found when the state is centered on a 
nodal line (soliton D for instance) or nodal point, but in 
this case it is unstable for low powers, developing lateral 
oscillations between two nodal positions (see Fig. 0Kf), 
simulations for different power D-type solitons). 

In conclusion, we have introduced and analyzed novel 
types of stable discrete solitons in two-dimensional pho- 
tonic lattices. Our simulations reveal that the exper- 
imental demonstration of our results is accessible with 
standard techniques. 
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